A flexible error estimate for the application of 
centre manifold theory 

Zhenquan Li A.J. Roberts* 

February 1, 2008 



Abstract 

In applications of centre manifold theory we need more flexible er- 
ror estimates than that provided by, for example, the Approximation 
Theorem 3 by Carr [Q, |6| . Here we extend the theory to cover the case 
where the order of approximation in parameters and that in dynam- 
ical variables may be completely different. This allows, for example, 
the effective evaluation of low-dimensional dynamical models at finite 
parameter values. 
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1 Introduction 



Interest in the dynamical behaviour of a physical system usually lies in the rel- 
atively low-dimensional evolution after heavily damped modes have become 
insignificant. There are many successful applications of centre manifold tech- 
niques to create models of these relatively simple dynamics. We here mention 



some applications to physical fluid mechanics. Iooss [|T6| , |T7f , and Laure (18 
analysed the dynamics of Taylor vortices in Taylor-Couette flow, whereas 
Chossat 0] and Hill discuss the non-axisymmetric dynamics involving 
mode competition by using centre manifold theory. Mode interactions in the 
dynamics of convection in porous media are analysed with centre manifolds 
by Neel |7|, || and Graham & Steen |TJ. Arneodo et al |], |, | re- 
duced the dynamics of triple convection down to a set of three coupled odes, 
numerically verified the modelling and then proved the existence of chaos. 
Roberts et al discussed centre manifolds of forced dynamical systems ||, and 
derived low-dimensional models using centre manifold techniques for contam- 
inant dispersion in channels |24| , shear dispersion in pipes [|5[] , thin film fluid 



dynamics [|33|] , coating flows over a curved substrate in space [^, 20], and 
Mei, Roberts & Li |22|, [23| derived models for turbulent shallow water flow 
written in terms of vertically averaged quantities derived from the k-e model 
for turbulent flow. Such applications assure us that centre manifold theory 
provides a useful route to the low-dimensional modelling of high-dimensional 
dynamical systems. 

Centre manifold theory guarantees the existence of low-dimensional mod- 
els, matches the solutions of original and the low- dimensional systems, and 
quantifies errors in the approximation. Algebraic techniques to construct 
low-dimensional models are based upon the theory. In problems specified 
in the standard form (D), the centre manifold may be calculated simply 
by iteration, see Carr [|J for example. For the more directly applicable 
form u = Cu + /(it, e), solutions in the form of an asymptotic power series 
are found using methods developed by Coullet & Spiegel || (and reinvented 
by Leen [19|]). The derivation of initial conditions for such low-dimensional 
models is given through projecting the initial condition of the system onto the 



centre manifold [|10] , 30] , [31, |34[]. But many of the applications require more 



flexible errors estimates. For example, physical models recoverd by evalua- 
tion at a finite value of a supposedly asymptotically small parameter often 
need high order approximations in the parameter PjL EI3, R3, EBL ETT 



Thus asymptotic errors estimates need to be made to high order in some 
parameters and only low order in other variables. In Section |3| we extend 
a theorem of Carr & Muncaster [HL |5|] to rigorously support such flexible 
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appr oximat ions . 

2 A simple example 

To introduce the issues, consider the well known prototype bifurcation prob- 
lem 

x = ex — xy , y = —y + x 2 , (1) 
where e is a parameter. By adjoining the trivial equation, it becomes: 

e = , x = ex — xy , y = —y + x 2 . (2) 

According to the distribution of the eigenvalues of the system, we may seek 
a centre manifold of form y = h(x, e). Substituting this into the system (|2]) 
we deduce that h must satisfy 

h = x 2 - ^-x(e - h) . (3) 
ox 

Solving this iteratively leads to the approximations 

fc(°)=0, h^=x 2 , =x 2 -2ex 2 + 2x\ etc. (4) 

Now elementary calculation shows that the above approximations hS""' sat- 
isfy (Q) to a residual O (|(e, x)\ n+2 ) as (e,x) — > 0; an error equivalently ex- 
pressed as O (e n+2 + x n+2 ) since a term ce p x q (for some constant c / 0) 
is O (e p + x q ) only if p'/p + q'/q > 1. Therefore the centre manifold is 
y = h(x, e) = + O (|(e, x)\ n+2 ) by, for example, Theorem 3 of g p264]. 

The limitation in applications is that the established theorem on approx- 
imation strongly couples the order of truncation in both parameters and 
variables — the "weight" of the parameter and variable is the same in the 
error estimate. Some flexibilty may be introduced by a nonlinear transfor- 
mation of the parameters; for example, introduce 5 = sfe and instead of (0) 
study 

5 = 0, x = 5 2 x — xy , y = —y + x 2 . (5) 

The resultant iterative solution of @ is identical to @. The only difference is 
that the approximation theorem asserts the errors in are O (\(5, x)\ 2n+2 ) 
as (5,x) — > 0, that is, O (e n+1 + x 2n+2 ). Thus certain trivial nonlinear trans- 
formations make no difference to the algebraic analysis but do affect the error 
estimate. There must be more flexibility in the errors than has so far been 
proved. 
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A more flexible error bound to include the effects of both e and x to 
any desired orders is to express and seek errors as 0(x q ,e p ). Note that 
/' — O (x q , e p ) means that any terms in / of the form ce p x q (for some constant 
c 7^ 0) must satisfy p' > p or q' > q. For example, we may deduce, supported 
by Theorem H herein, 



h = x 2 (l - 2e + Ae z ) + x\2 - 16e + 88e") + O [x\ e 6 ) . (6) 

This kind of error allows us separately to choose the orders of the parameter e 
and the variable x which we want to include in the centre manifold. For 
example, here we may compute to higher orders in e, observe the pattern 
of coefficients in this simple problem, and realise |29| that the above is the 
low-order Taylor polynomial in e of 

Such approimations to high-order in parameters and low-order in dynamical 
variables were used, before proof, and are essential to the analyses in [J^, 
E2l BJS, E3LE1. 



3 The flexible extension 

Consider dynamical systems expressed in the form 

x = Ax + f(x,y,e) , 

y = By + g(x,y,e) , (8) 

where x e R m , y £ R n , e £ R/ and A and B are constant matrices such 
that all the eigenvalues of A have zero real parts, and all eigenvalues of 
B have negative real parts. Functions / and g are nonlinear for x, y, e 
and /(0,0,0) = g(0,0,0) = /'(0,0,0) = <?'(0,0,0) = (where /' = 
[fx, fyi fe\i an d similarly for g' and other Jacobians). 

For any function cf) : R m xR'-> R n which is a continuously differentiable 
function and </>(0, 0) = <p'(0, 0) = 0, define 

(H<f>) = (f> x ( x , e)[Ax + f(x, 4>{x, e), e)] - B<j>(x, e) - g(x, <f>(x, e), e) . 

Also let O (s g , e p ) denote any terms of the form cef 1 . . . ef'sf 1 . . . s^ 1 (where 
the constant c ^ 0) which satisfy p\ + • • • + p\ > p or q x + • • • + q m > q and 
PiAj > 0. 
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Theorem 1 (Approximation) Suppose that 

(H<j>) (x, e) = 0(x q , e p ) as (x, e) -»■ , 

where p > I, q > 1, then 

\h(x,e) - <f>(x,e)\ = 0(x q ,e p ) as {x,e)^0. 

That is, the errors in the approximation <fr to a centre manifold h is the same 
as the order of the residuals of the equations of the dynamical system. 

Proof: This proof is adapted from Carr ||, pp25-28]. 

Let 6 : R m x K l — > R n be a continuously differentiable function with 
compact support such that 



Set 



6(x, e) = (f>(x, e) for |(a;, e)| small. 



N(x,e) = x {x, e) [Ax + F (x, 6(x, e), e)] 

-BG(x,e) -G(x,6(x,e),e) , (9) 



where 



F(x,y,e) = f (xip(^j ,y,e 

G(x,y,e) = g(x^(jj,y,e 

where ip : R m — > [0, 1] is a infinitely differentiable function with vjj(x) = l 
when | a; | < 1 and ip(x) = when \x\ > 2 and S is a positive real number. 
The properties of F and G are the same as in j|, pl8] for e small. So 
JV(aj, e) = O (x 9 , e p ) as (aj, e) -> 0. 

For a > and 6 > let T be the set of Lipschitz functions : R m x 
R' — >• R n with Lipschitz constant b, \h(x,e)\ < a for (x, e) G R m x R' and 
/i(0, 0) = 0. With the supremum norm || . ||, T is a complete space. 

For ZieT and a?o G R m , let x(t,Xo,h) be the solution of 

x = Ax + F(x, h, e) , x(0, x , h) = x . 

The bounds on F and h ensure that the solutions of the above equation 
exists for all time t. Define an operator T on T by 

(Th)(x ) = / e~ S G (x(s, x , h), h(x(s, x , h)), e) ds . 
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We know that T is a contraction mapping and the centre manifold y = 
h(x, e) is a fixed point of T for a, b and 5 small enough from the proof of 
the existence theorem [HI pp. 16-19]. Define 

SZ = T(Z + 6) , 

for Z such that Z + 6 G T. The domain of S is a closed subset of T since 
G r. Since \SZ X - SZ 2 \ = \T(Z 1 + 6) - T(Z 2 + 6)1 thus S is also a 
contraction mapping. For K > letQ 

y = |Z G T | |Z(x,e)| < K(x q ,(?),V(x,e) gOcCxR 1 }, 

where O is a neighbourhood of the origin in R m x R'. Since iV(a;,e) = 
0(x q , e p ) as (x, e) -> 0, then 

|JV(aj,e)| < Ci(x 9 ,e p ), (x,e)eO (10) 

where Ci is a constant. Thus Y is not empty because N G y C T by defining 
0(a;, e) such that Ci < if. If we can find a constant K such that 5 maps Y 
into y, then 3 Z G F is a fixed point of 5*, and 

Zo=S(Z )=T(Z + 0)-0, 
that is, T(Z + 0) = Z + 6 , 

i.e., -Zo + 6 is a centre manifold of (Bp, let h = Zq + 6, 

\h(x,e) - 0(x,e) \ = Z < K(x\e p ) . 

To finish the proof define 

Q(x,Z,e) = e x (x,e)[F(x,e + Z,e)-F(x,e,e)]-N(x,e) 
+ G(x,0 + Z,e) - G(x,0,e) . 

Then 

\Q(x,Z,e)\ < \Q(x,0,e)\ + \Q(x,Z,e) -Q(x,0,e)\ 

= \N(x,e)\ + \Q(x,Z,e) - Q(x,0,e)\ . (11) 

From the properties of F and G on pl8 in [HJ and #'((), 0) = 0, we have 

\Q(x,Z,e)-Q(x,0,e)\<k(5)\Z\ for \Z\ < 5 . (12) 

1 K(x q ,e p ) denotes the product of K and sum of finite terms ce^ 1 ■■■( p l l x q ^ ■ ■ • , 
where pi + • • • + Pi > V or 1\ + ' ' ' + ?m > 5 an d Pi,Qj > 0, c is a non-zero constant. 
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Using (0), Q and (0), 

\Q(x,Z,e)\ < (Ci + A'£;(5))(x 9 ,e p ), for Ze7. (13) 

Using the same calculations as (2.5.9) on p27 M, for each r > 0, there is a 
constant M(r) such that 

\x(t,x ,e)\ < M(r)\x \e-^ , t<0 (14) 

where 7 = r + 2M(r)/c(<5) and x(t, x , e) is the solution of 

x = Ax + F(x, Z(x, e) + 0(x, e), e) , x(0, a? , e) = a; . 

Using (2.3.6) on pl8 and (2.5.3) on p26 in |§, and ([Lj), ([14]), if Z G F 



|(5Z)(a;o,e) < 

< 



e- Bs (d + Kk(5))(x q ,e p ) ds 

-00 




e-^(Ci + ia(5))M(r)"e- 97S (4, e p ) 

< C(d + If fc(*))M(r)9(/3 - 7g) _1 (4, e P ) • 

provided 5 and r small enough so that (3 — > 0. Choose O and 5 small 
enough and K large enough such that K > C(C\ + Kk{8))M{r) q ((3 — 7<?)~ 1 . 
Therefore S : Y — >• F . Hence Theorem |I] holds. 4 

More general theorems allowing varying orders of truncations within pa- 
rameters and dynamical variables may be also useful. However, most such 
cases can be easily established by simple nonlinear transformations of the 
parameters along the same lines as the example in (^). 

In applications, such as many fluid dynamics problems, we need the- 
ory not only dealing with infinite dimensional problems, but also infinite 
dimensional centre manifolds. Carr presented the corresponding results 
for infinite dimensional problems, and analysed two problems arising from 
partial differential equations for finite dimensional centre manifolds. The 
restrictions upon the nonlinear terms / is that / has order 2 continuous 
derivative and /(0) = /'(0) = 0. More recently, Gallay gave an exten- 



sion of the existence theorem to infinite dimensional centre manifolds, but a 
bounded restriction on the nonlinear terms is required. This condition limits 
its rigorous application. Scarpellini |36| apparently places significantly less 



restrictions upon the nonlinearities in the dynamical equations, but while he 
addresses infinite dimensional centre manifolds, the results are severely con- 
strained by requiring finite dimensional stable dynamics. Haxagu§ [|14], [HJ 
has developed theory supporting infinite dimensional models, such as the 
Korteweg-de Vries equation, but only by placing extreme restrictions upon 
the linear operators. We identify the extension of the theorems to infinite 
dimensional centre manifolds as a significant problem for future research. 



AJ Roberts, February 1, 2008 



References 



8 



Acknowledgement: we thank the University of Southern Queensland and 
the Australian Research Council for supporting this research. 

References 

[1] A. Arneodo, P.H. Coullet, and EA. Spiegel. The dynamics of triple 
convection. Geophys. Astro. Fluid Dyn., 31:1-48, 1985. 

[2] A. Arneodo, P.H. Coullet, EA. Spiegel, and C. Tresser. Asymptotic 
chaos. Physica B, 14:327-347, 1985. 

[3] A. Arneodo and O. Thual. Direct numerical simulations of a triple 
convection problem versus normal form predictions. Physic Letters A, 
109:367-372, 1985. 

[4] J. Carr. Applications of centre manifold theory, volume 35 of Applied 
Math Sci. Springer- Verlag, 1981. 

[5] J. Carr and R.G. Muncaster. The application of centre manifold the- 
ory to amplitude expansions. I. ordinary differential equations. J. Biff. 
Eqns., 50:260-279, 1983. 

[6] J. Carr and R.G. Muncaster. The application of centre manifold theory 
to amplitude expansions. II. infinite dimensional problems. J. Biff. Eqns, 
50:280-288, 1983. 

[7] P. Chossat, Y. Demay, and G. Iooss. Competition of azimuthal modes 
and quasi-periodic flows in the Couette- Taylor problem. In Bynamics of 
infinite- dimensional systems, volume 37 of NATO Adv. Sci. Inst. Ser. 
F: Comput. Systems Sci., pages 45-56. Springer, 1987. 

[8] P.H. Coullet and EA. Spiegel. Amplitude equations for systems with 
competing instabilities. SIAM J. Appl. Math., 43:776-821, 1983. 

[9] S.M. Cox and A.J. Roberts. Centre manifolds of forced dynamical sys- 
tems. J. Austral. Math. Soc. B, 32:401-436, 1991. 

[10] S.M. Cox and A.J. Roberts. Initial conditions for models of dynamical 
systems. Physica B, 85:126-141, 1995. 

[11] Th. Gallay. A center-stable manifold theorem for differential equations 
in banach spaces. Commun. Math. Phys, 152:249-268, 1993. 



AJ Roberts, February 1, 2008 



References 



9 



[12] M.D. Graham and P.H. Steen. Strongly interacting traveling waves 
and quasi-periodic dynamics in porous medium convection. Physica D, 
54:331-350, 1992. 

[13] A. Hill and I. Stewart. 3-mode interactions with o(2) symmetry and a 
model for Taylor- Couette flow. Dynamics Stability Systems, 6:267-339, 
1991. 

[14] M. Haragu§. Model equations for water waves in the presence of surface 
tension. Eur. J. Mech, 15:471-492, 1996. 

[15] M. Haragu§. Reduction of pdes on unbounded domains application: 
unsteady water waves problem. J. Nonlinear Sci., 8:353-374, 1998. 

[16] G. Iooss and M. Adelmeyer. Topics in Bifurcation Theory. World Sci, 
1992. 

[17] G. Iooss and J. Los. Bifurcation of spatially quasi-periodic solutions in 
hydrodynamic stability problems. Nonlinearity, 3:851-871., 1990. 

[18] P. Laure and Y. Demay. Symbolic computation and equation on the 
centre manifold: Application to the couette-taylor problem. Computers 
And Fluids, 16:229-238, 1988. 

[19] T.K. Leen. A coordinate independent center manifold reduction. Phys. 
Letts A, 174:89-93, 1993. 

[20] Zhenquan Li and A.J. Roberts. The accurate modelling of thin 3D fluid 
flows with inertia on curved substrates. In Proceedings of EMAC98, 
pages 315-318, 1998. 

[21] Zhenquan Li and A.J. Roberts. The accurate and comprehensive model 
of thin fluid flows with inertia on curved substrates. Technical report, 
| http : //xxx . lanl . gov/abs/chao-dyn/9906011|| , June 1999. submit- 
ted to IMA J Appl Math. 

[22] Z. Mei and A.J. Roberts. Equations for turbulent flood waves. In 
A. Mielke and K. Kirchgassner, editors, Structure and dynamics of non- 
linear waves in fluids, pages 342-352. World Sci, 1995. 

[23] Z. Mei, A.J. Roberts, and Zhenquan Li. Modelling the dynamics 
of turbulent floods. Technical report, [ tittp : / /xxx . lanl . gov/abs/ 



patt-sol/9907001|| , June 1999. submitted to SIAM J Appl Math. 



AJ Roberts, February 1, 2008 



References 



10 



[24] G.N. Mercer and A.J. Roberts. A centre manifold description of con- 
taminant dispersion in channels with varying flow properties. SIAM J. 
Appl. Math., 50:1547-1565, 1990. 

[25] G.N. Mercer and A.J. Roberts. A complete model of shear dispersion 
in pipes. Jap. J. Indust. Appl. Math., 11:499-521, 1994. 

[26] M.C. Neel. Convection in a horizontal porous layer of infinite extent. 
European J. Mech. B Fluids, 9:155-176, 1990. 

[27] M.C. Neel. Natural convection in a horizontal porous layer of infinite 
extent: inhomogeneous heating. In European Conference on Iteration 
Theory, pages 272-278. World Sci. Publishing, 1991. 

[28] M.C. Neel. Nonhomogeneous data in the choice of the convective con- 
figurations in a porous layer. Nonlinear Vib. Probl., 25:297-311, 1993. 

[29] A.J. Roberts. Simple examples of the derivation of amplitude equations 
for systems of equations possessing bifurcations. J. Austral. Math. Soc. 
B, 27:48-65, 1985. 

[30] A.J. Roberts. The utility of an invariant manifold description of the 
evolution of a dynamical system. SIAM J. Math. Anal., 20:1447-1458, 
1989. 

[31] A.J. Roberts. A sub-centre manifold description of the evolution and in- 
teraction of nonlinear dispersive waves. In L. Debnath, editor, Nonlinear 
waves, chapter 9, pages 127-156. World Sci, 1992. 

[32] A.J. Roberts. Low- dimensional models of thin film fluid dynamics. Phys. 
Letts. A, 212:63-72, 1996. 

[33] A.J. Roberts. An accurate model of thin 2d fluid flows with inertia on 
curved surfaces. In P. A. Tyvand, editor, Free- surface flows with viscosity, 
volume 16 of Advances in Fluid Mechanics Series, chapter 3, pages 69- 
88. Comput Mech Pub, 1998. 

[34] A.J. Roberts. Computer algebra derives correct initial conditions for 
low-dimensional dynamical models. Comput. Phys. Comm., 1999. to 
appear. 

[35] R. Valery Roy, A.J. Roberts, and M.E. Simpson. A lubrication model of 
coating flows over a curved substrate in space. Technical report, [http : 



//xxx. lanl .gov/abs/patt-sol/9705002|| , August 1997. submitted to 



J Fluid Mech. 



AJ Roberts, February 1, 2008 



References 11 

[36] B. Scarpellini. Centre manifolds of infinite dimensionals. I. Main results 
and applications. Z. Angew. Math. Phys., 42:1-32, 1991. 



AJ Roberts, February 1, 2008 



